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In contrast to the experimentally widely used isentropic de- 
magnetization process for cooling to ultra-low temperatures 
we examine a particular classical model system that does not 
cool, but rather heats up with isentropic demagnetization. 
This system consists of several magnetite particles in a col- 
loidal suspension, and shows the uncommon behavior of disor- 
dering structurally while ordering magnetically in an increas- 
ing magnetic field. For a six-particle system, we report an 
uncommon structural transition from a ring to a chain as a 
function of magnetic field and temperature. 

05.70.Ce, 05.70.Fh, 05.70.Jk, 75., 75.10., 75.50. 



Cooling to ultra-low temperatures is presently exper- 
imentally achieved using the isentropic demagnetization 
process suggested in 1926 by Debye and Giauque [|l],D. 
Here we study the thermal and magnetic properties of 
a particular model system consisting of six magnetite 
nanoparticles in a colloidal suspension. Our calculations 
indicate that changing the temperature or magnetic field 
leads to a transition from a ring-like to a chain-like struc- 
ture. We demonstrate that such a model system, that 
could be realized in a ferrofluid, would not cool but rather 
heat up during an isentropic demagnetization. The sys- 
tem of magnetic nanoparticles shows the uncommon be- 
havior of disordering structurally, i.e. increasing its vol- 
ume in phase space, while ordering magnetically in an 
increasing magnetic field. This behavior, associated with 
the break-up of the relatively rigid rings to floppy chains, 
occurs once the energy gain upon aligning a chain of 
dipoles with the field exceeds the energy cost of break- 
ing up a ring. Thus systems that disorder structurally in 
high fields may be used as coolants based on isentropic 
magnetization instead of demagnetization. So far, only 
ferrofiuid systems with a large number of particles have 
been discussed as candidates for use in magneto-caloric 
heat engines H. Here we report a paradoxical magnetic 
cooling phenomenon that is unique to systems with only 
few particles. 

Both the conventional and paradoxical process utilize 
the energy change associated with particular structural 
changes for cooling. The conventional isentropic demag- 
netization process uses the fact that a magnetic system, 
such as a spin system, orders magnetically and thus low- 
ers its entropy in presence of an external magnetic field. 



Removal of the external magnetic field at constant tem- 
perature leads to an increase in entropy due to magnetic 
disordering, which requires energy. Decreasing the exter- 
nal magnetic field at constant entropy consequently leads 
to a temperature lowering. With this method, systems 
such as copper have been cooled down to temperatures 
as low as 50 nK Q . 

In the following we address model systems consisting 
of a finite number of magnetic particles that may undergo 
structural transitions. Both the structural and magnetic 
degrees of freedom are important in this system and can 
not be decoupled. Applying a sufficiently high magnetic 
field causes the system to order magnetically while disor- 
dering structurally, at the cost of internal energy. Conse- 
quently, such a system will exhibit the paradoxical phe- 
nomenon of cooling by isentropic magnetization. 

The transformation from a ring to a chain is associ- 
ated with freeing up a structural degree of freedom, with 
a corresponding increase in entropy. Such a transfor- 
mation can be induced by a magnetic field in a system 
of magnetic dipoles, where the energetics is governed by 
dipole-dipole interactions between the particles and an 
interaction of each particle with the external field. In 
small fields, the ring is stabilized with respect to the 
chain structure if the gain in dipole-dipole interaction 
upon connecting the chain ends energetically outweighs 
the dipole misalignment energy in a bent structure. The 
energy gain upon aligning all individual dipoles with a 
sufficiently high applied field, on the other hand, stabi- 
lizes the chain structure. 

One system that satisfies all the requirements on a 
paradoxical magnetic coolant consists of a few (4<iV< 14) 
super-paramagnetic particles of magnetite. Such parti- 
cles are the key constituents of ferrofiuids, which have 
attained rapidly increasing interest during the past few 
years |^,|[. Recently we have shown that such systems 
exhibit intriguing phase transitions between the ordered 
ring and chain phases and one disordered phase [Q. We 
also pointed out that self-assembly in these systems could 
be used to store information [p|. 

In the following, we study the thermodynamic behav- 
ior of a six-particle system, where the chain and the ring 
are the dominant stable structural isomers. We chose 
this particular system as it allows a simple discussion of 
the two thermodynamical features characteristic of iso- 
mer transitions in finite systems, namely the isomeric 



phase space and the transition probabihty that is hnked 
to the transition time. It is true that such ring and chain 
isomers are also stable at much larger sizes. Nevertheless, 
the thermodynamical behavior becomes more complex in 
larger systems, where the number of relevant structural 
isomers grows rapidly with increasing N. 

Here we calculate the entropy and the dependence of 
temperature on an external magnetic field at constant 
entropy {dT/dB)s for a model system of six magnetite 
particles with a radius cr = 50 A and a permanent mag- 
netic moment oi ji — 2.63x10^ fx-Q. The total potential 
energy U of this system is given by l6l 

N 
U = Y1 {(^o/4) [A* • Ai - 3(A« • ^y)(Aj • f^J)] (1) 
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The pairwise interaction energy is given by the dipole- 
dipole and a non-magnetic interaction energy. The latter 
is dominated by a repulsion between the spherical par- 
ticles, but contains also a weak attractive part due the 
van der Waals interaction. It is modelled by the above 
Morse- type potential with parameters [Q e = 15.1 /zeV 
and p = 2.5 A. The second sum reflects the interaction 
between the magnetite particles and the external mag- 
netic field -Boxt that is aligned with the z axis. 

We would like to point out that the finite ferrofluid 
systems undergo all the intriguing transitions described 
below occur independent of these parameter values. Our 
particular choice has been taken to bring the transition 
into an experimentally accessible and interesting region 

All thermodynamic quantities can be derived from the 
canonical partition function Z(B,T) by appropriate dif- 
ferentiation. We determined Z using the Metropolis 
Monte Carlo method ||ll|, which we combined with a 
special type of optimized data analysis p2| to calculate 
all thermodynamic properties as functions of the temper- 
ature T and the magnetic field i?oxt as external variables 
. The entropy S is given by 



S = kBi HZ) - P 



d\n{Z) \ 



(2) 



where f3 — l/fceT. The fundamental thermodynamic 
expression dE = TdS — ^zdB yields immediately the 
Maxwell relation 



dT 



- {dS/dB)j, 
{dS/dT)g 



(3) 



that describes the temperature response to the external 
field in isentropic processes. We calculate this quantity 
using the expectation values of the potential energy U 
and the z-component of the magnetic moment /i^ of the 
whole system as 



dT 
dB 






(4) 



Our results for S and (§5)0 as a function of Soxt and 
T are represented by contour plots in Fig. yj and Fig. y, 
respectively pSl]. 
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FIG. 1. Contour plot of the entropy S' as a function of the 
external magnetic field Bcxt and temperature T. The model 
system discussed here shows a temperature decrease of up 
to few degrees Kelvin as the field B is increased at constant 
entropy. 

The steps in the isentropes displayed in Fig. ^ indi- 
cate a temperature decrease with increasing magnetic 
field at constant entropy. For example, an increase of 
the field from 45 to 55 Gauss cools the system down from 
T = 12 K to T^8 K. The narrow region in the Boxt - T 
space, where these kinks occur, separates the chain and 
ring phases. It is only in this narrow region of the B^y^t—T 
space that {dT/dB)s shows a nonzero value and hence 
a potential for magnetic cooling, indicated in Fig. g. A 
closer inspection of Fig. |l| shows that the S = const, 
lines change their slope in the chain phase at high fields, 
which also corresponds to a positive value of {dT/dB)s 
in Fig. H This behavior is simply related to the fact 
that chains behave like a conventional magnetic system, 
since their magnetic and structural order increases with 
increasing external magnetic field. In the six-particle sys- 
tem discussed above, the conventional cooling mechanism 
by isentropic demagnetization in high fields is about one 
order of magnitude less important than the paradoxical 
magnetic cooling that is related to the ring-chain transi- 
tion. 

Let us now discuss some details related to a possible 
realization of the system. In our simulations we disre- 
garded the internal degrees of freedom of the magnetic 
particles and the colloidal suspension. Obviously, the 
cooling efficiency of this particular composite system is 
highest when the magnetic degrees of freedom of the par- 
ticles dominate. A potential candidate paradoxical cool- 
ing system is a dilute gas consisting of both magnetic 
particles, which could aggregate, and non-magnetic par- 
ticles of similar size. The actual thermal transition is ini- 



tiated by vibrational excitations that lead either to the 
breakup of rings into chains or a reconnection of chain 
ends to a ring. 
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FIG. 2. Contour plot of {dT/dB)s as a function the ex- 
ternal magnetic field Bext and temperature T. 

Let us now consider a mixture of chains, rings, and 
other isomers exposed to an oscillating magnetic field. 
Whereas the breakup energy increases as the neighbor- 
ing dipoles gradually lign up in the larger rings, the vi- 
brational excitation spectrum gets dense especially in the 
range of low frequencies that drop with I/7V^. Hence it 
is preferentially the larger rings that are broken up to 
chains by a low-frequency external field. Also the chain 
structures experience an analogous softening of the vi- 
brational spectrum with the factor of l/N"^. Yet due to 
the larger phase space volume, the likelihood of two chain 
ends to reconnect to a ring decreases drastically in larger 
systems. 

The chain ends (which are to meet for a chain-to-ring 
transition) not only move more slowly with larger N but 
have to explore a much larger space in which the other 
end can be. The time required for the two ends to re- 
connect rises dramatically with N , so that chain-to- ring 
transitions can not be realized in a reasonable time in an 
ensemble of long chains. 

In a realistic experiment, we also have to consider that 
the average size of the ferrofluid aggregates depends not 
only on the initial conditions, but increases as a func- 
tion of time due to further aggregation. To prevent this 
gradual size increase, we propose to expose the ferrofluid 
suspension to a low frequency magnetic field that should 
preferentially break apart the larger aggregates and to 
stabilize the mixture of prevalent isomers in the region 
of five to seven particles per aggregate. 

Paradoxical magnetic cooling is by no means restricted 
to the model system discussed here. We believe that the 
same effect should also occur in other nanostructures, 
such as transition metal clusters, and even in bulk matter 
consisting of finite-size substructures. We hope that our 
results will stimulate a search for experimentally realiz- 
able systems that may find application in the fascinating 
field of ultra-low temperature physics. 
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